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Abstract 



Propagation of initially localized perturbations is investigated in chaotic cou- 
pled map lattices with long-range couplings decaying as a power of the dis- 
tance. The initial perturbation propagates exponentially fast along the lattice, 
with a rate given by the ratio of the maximal Lyapunov exponent and the 
power of the coupling. A complementary description in terms of a suitable 
comoving Lyapunov exponent is also given. 
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Propagation of fronts in spatially extended systems is a topic of wide interest, in several 
scientific contexts including fluid dynamics, dendritic growth, directional solidification, liquid 
crystals, chemical reactions, flame propagation and epidemics spreading In one spatial 

dimension, partial differential equations 0, coupled map lattices (CML) p and cellular 
automata [^] models have been studied. The fronts were either separating stable and unstable 
steady states or regular and chaotic regions ||lO| , p |. 



A similar phenomenon is the propagation of disturbances in fully chaotic states of CMLs 



T^ , p!3[| and of the complex Ginzburg-Landau equation |jT^. When the spatial coupling is 



local, e.g. for CML with nearest-neighbour interactions, two distinct regimes have been 



found to characterize the dynamics of fronts |T2|. In the first regime, the velocity of the 
propagating front is calculated in the framework of a linear analysis, while in the second the 
value of the velocity is determinated by the full nonlinear evolution of the system. For a 
particular class of models, a transition between the two regimes has been observed when a 
parameter is varied . The scenario is not expected to be different for couplings extending 



to more than two neighbours, since this affects only the limit velocity of propagation [|T5|. 

If more general spatial couplings are considered, even the very definition of a velocity can 
be nontrivial. For example, epidemic models in one dimension exhibit a finite propagation 
speed if the spatial coupling (i.e. the infection rate) decays exponentially or faster [Q]. The 
same result for some CML models is reported in Ref . [|16| , where it is observed that the time 
for a localized disturbance to overcome a threshold value at a certain distance I from the 
initially perturbed site grows linearly with / if the coupling is local or exponentially decaying, 
while is almost independent of / for couplings decaying with power laws (with exponent not 
too large), indicating that the velocity is "infinite". 

This Letter focuses on the study of disturbance propagation in systems with long range 
coupling, whose strength decays as a power law in space. The spatio-temporal evolution of 
an initially localized perturbation of a chaotic state is studied theoretically and numerically. 
The perturbation is found to spread exponentially fast along the lattice and an expression for 
the corresponding rate is given. A comoving Lyapunov analysis confirms that the prediction 
is indeed correct. 

In order to mimic a spatially extended system with long range interaction, we introduce 
the CML model 



uix,t + l) = f{u{x,t)) (1) 

where the indices x labels the sites of a chain of length L and t is the discrete time variable. 
The function f{x) is a chaotic map of the interval and 

u{x,t) = {1- e)u{x,t) + e J2 w{x-y)u{y,t) , (2) 

f)<\x-y\</^L 

where = {L is assumed to be odd). As usual, periodic boundary conditions are 
considered u{x,t) = u{x ± L,t). The coupling constant e is bounded between and 1, and 
we consider the coupling strength to decay in space as 

/ \ c(a) 

w(x) = (3) 
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where a > 1 to insure that the normahzation constant c(a) = [J2y is bounded. 

Obviously, this constant is independent of L in the thermodynamic hmit L ^ oo. For 
a — i> +00 the model (El) reduces to the usual nearest-neighbours CML p[. Notice that, 



at variance with the model of Ref. [|16[, the long range coupling is not introduced as a 
perturbation of the nearest-neighbours CML. 

To study the propagation of localized disturbances in the system (H), let us consider 
two chaotic trajectores {u{x,t)} and {v{x,t)} generated by starting initial conditions which 
differ only around a single site x = 0. More precisely, we assume that u{x, 0) is a typical 
chaotic state, obtained after all transients have died out, and v{x, 0) = u{x, 0) +Uo{x), where 
Uq{x) is a function localized around the origin. 

If only linear mechanisms of propagation are present |]T2| , |T3| , we can assume that the 
evolution of 6u{x,t) = u{x,t) — v{x,t) is ruled by the linearized dynamics 



6u{x, t 



m{x, t) 



[1 — e)6u{x,t) + e ^ w{x — y)6u{y,t) 



0<\x-y\<AL 



(4) 



where m{x, t) = f'{u{x, t)) is the local multiplier along the assigned trajectory. This hypote- 
sis is justified at least for large distances ^ 1) where the disturbance is small. Moreover, 
its validity has also been numerically checked by comparing the evolution given by Eq. 
with that of u{x,t) — v{x,t). As a matter of fact, nonlinear effects, i.e. saturation of the 
perturbation growth, are present only in the central part of the disturbance. 

We consider now the spatial shape of the leading edge of the front and its temporal 
evolution in the limit t ^ 00 for the ideal case of an infinite lattice [L = 00). As a starting 
point, we make the ansafe that for ^ 1 



6u{x, t) 



\x\ 



(5) 



with (t){x,t) weakly dependent on x. Inserting Eq. (|^) in Eq. (^), we obtain 



+ 1) 



\x 



1/3 



m(x, t) 



\x\ 



(6) 



Due to the long range coupling we can assume that, as a first approximation, a mean field 
description holds. This amounts to neglect the spatial fluctuations of t) and to replace 
m{x,t) with its average e^, where A is the (maximal) Lyapunov exponent. Eq. (^ is then 
splitted in two equations, one for the time evolution 



0(t + l) =eV(t) 



and one for the spatial profile 



\x\ 



E 

yj^x 



W[X 



y) 



\y\ 



(7) 



(8) 



Moreover, at least for x ^ 1, the sum appearing in Eq. (^ can be approximated by neglecting 
the small terms coming from decaying tails, i.e. by extending the sum only between 1 and 
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X — 1 (for simmetry reasons, we can also consider x > 0). The spatial shape of the front is 
thus conserved in time if the fixed-point condition 



1~ V ^(^) (g) 

is satisfied for x large enough. Since the leading contributions to the sum came from the 
extrema and can be estimated to be of order \x\~^ and |x|~°, Eq. (|^) is fulfilled for /5 < a. 
By combining this result with that of Eq. (0), we get the expression for the asymptotic 
behaviour of the front leading edge 

6uix,t)r.'-^ . (10) 

Let us define the front position r{t) as the maximal distance from a; = where 
\Su{x,t)\ > 9, with 9 > being a preassigned threshold. According to Eq. (|1Q|) , r{t) 
grows exponentially as 



r{t) ~ exp -t = exp {S{P)t) (11) 




for t — >■ oo, where the value of the (3 parameter depends on the initial shape of the dis- 
turbance. From the above reported arguments it is clear that, if the perturbation decays 
initially as Uq{x) ~ with (3 < a, than its profile will be preserved during the time 

evolution. For more general initial conditions the selected (3 can be determined by means 
of the following argument. Let us consider uo{x) to be a superposition of several profiles, 
each one decaying as but with different (3 < a. In the linear approximation, each one 

of these components will propagate independently with a different rate S{f3). On the basis 
of general arguments we expect that (for Uo{x) sufficiently localized) the profile with 
slowest growth rate S{a) = A/a will be selected. 

These results are analogous to those found for CML with nearest-neighbours coupling 
p!3| and for the complex Ginzburg-Landau equation the main difference being that the 
tails of the front are exponential for the latter case. 

We numerically tested the above predictions for a lattice of coupled tent maps (/(-z) = 1 — 
2\z\), with several values of the coupling constant and chain lengths ranging from L = 1001 to 
L = 20001. In particular, we have considered a single-site perturbation uq{x) = 6^^, where 
y is the usual Kroneker delta. We computed the time evolution of 6u{x,t) according 
to Eq. (^ and averaged it over different reference trajectories. As shown in Fig. 1, the 
perturbation profile decays on average as while growing exponentially in time. We 

have also verified that an initial disturbance with a decaying profile \x\~^ is conserved only 
for (3 < a, otherwise a power law decay exponent a is always found. Finally, we checked the 
validity of Eq. (|^) by verifying that, for increasing t, the local growth rate on the tails of Su 
approaches the Lyapunov exponent. 

A direct numerical test of Eq. (|TT|) is complicated by the presence of finite-size effects 
(see however Ref. for results on a closely related CML model). An indirect check is 



accomplished by considering a suitably defined comoving Lyapunov exponent A |T^. For 
CML with local coupling, A{v) is defined as the asymptotic growth rate of a disturbance in 
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a reference frame moving along the "world line" Xy{t) = vt, where v is the frame velocity 
(1^1 < 1). This amounts to assume Su{xy(t),t) ~ exp(A(f)t). Within this scheme, the 
condition A{v) = defines the propagation speed of an initially localized disturbance |T8| , p!3 



In the present case, since the long range coupling leads to exponentially fast propagation 
(see Eq. (0)), we rather define comoving Lyapunov exponent A{R) in a reference frame 
moving along the "world line" XR{t) = [exp(i?t) — 1]. Therefore, on the leading edge and 
for sufficiently long time t, the relation 

6u{x, t) ~ exp (A(_R)t) = exp (At — a log x) (12) 

should hold. From this equation it is readily seen that the comoving exponent must be a 
linear function of the rate namely 

A(i?) = \-aR . (13) 

Notice that, in analogy with system with local coupling, the condition A{R) = gives 
exactly the growth rate S{a) = A/a, in agreement with the above prediction. 

As before, the finite size of the system prevents the numerical computation of A for 
asymptotically large times as requested by its very definition. Therefore, we computed its 
finite-time value 



A(i?,t) = -(log 



6u{xR(t),t) 



6u{0,0) 



(14) 



where R = log |x + and (■) is the average over different reference trajectories. As can be 
easily realized, the maximal accessible rate will decrease as log L/t. Therefore if the iteration 
time is doubled, the chain length should be increased by a factor L in order to achieve the 
same maximal R. Due to the large of amount of CPU time required by the iteration, it is 
thus not feasible to consider system sizes larger than 10^, and the accessible ranges of R and 
t values are limited by this constraint. Nevertheless, the results reported in Fig. 2 confirm 
that a linear behaviour A{R,t) = X*{t) — aR is observed at each time and for R not too 
small, but with an intercept A* 7^ A due to the finite time of observation. Empirically, we 
found that this intercept converges according to the rule X*{2nt) — A*(t) ^ constant/g""^, 
for some q > 1, so that we can extrapolate its asymptotic value X*{oo) on the basis of the 
available data. Indeed, for the situation reported in Fig. 2 (a = 3/2 and e = 1/3) we have 
estimated a value A* (00) = 0.335, in excellent agreement with the corresponding Lyapunov 
exponent A = 0.338. 

The deviation from Eq. (|13D at small values of R is due to transient effects. We have 
numerically observed that, for increasing t, the interval of R values where deviations are 
observed reduces. 

In conclusion, we have fully identified the mechanism that rules the disturbance propa- 
gation for systems with power law long range couplings. This gives at any position in space 
an exponential increase in time, and a power law fall-off with x. The power with which 
the perturbation decays, is, for generic initial conditions equal to the power describing the 
interaction fall-off. Moreover, the time needed for the disturbance to propagate with finite 
amplitude at a given distance / is inversely proportional to the Lyapunov exponents and 
increases logarithmically with I. 
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FIGURES 



FIG. 1. Plot of the logarithm of the average disturbance amplitude {\Su{x,t)\) versus Inx at 
three different times for coupled tent maps with a = 2, L = 5001, e = 1/3. From bottom to top 
the three solid curves correspond to t = 20,40 and 60, respectively. The dashed line has a slope 2. 

FIG. 2. The finite-time comoving Lyapunov exponents A(i?, t) as a function of rate R at time 
t = 20,40 and 80 (from bottom to top), for coupled tent maps with a = 3/2, L = 5001, £ = 1/3. 
The dashed line represents the asymptotic expression A(i?) = X — 3/2 ■ R. 
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